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ABSTRACT 



Scaling violation of inclusive jet production at small- X in hadron scattering with increasing total 
collision energy is discussed. Perturbative QCD based on the factorisation theorem for hard 
processes and GLAPD evolution equations predicts a minimum for scaled cross-section ratio 
that depends on jet rapidity. Studies of such a scaling violation can reveal a vivid indication of 
new dynamical effects in the high-energy limit of QCD. The BFKL effects, which seem to be 
seen in recent L3 data at CERN LEP2, should give different results from GLAPD predictions. 
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QCD is an essential ingredient of the Standard Model, and it is well tested in hard processes 
when transferred momentum is of the order of the total collision energy (Bjorken limit: Q 2 ~ 
s — > oo). The cornerstones of perturbative QCD at this kinematic regime (QCD-improved 
parton model): factorization of inclusive hard processes |IJ and the Gribov-Lipatov-Altarelli- 
Parisi-Dokshitzer (GLAPD) evolution equation 0] provides a basis for the successful QCD- 
improved parton model. The factorisation theorem Q for inclusive hard processes ensures 
that the inclusive cross section factorises into partonic subprocess(es) and parton distribution 
function(s). The GLAPD evolution equation governs the log Q 2 -dependence (at Q 2 — > oo) of 
the inclusive hard process cross-sections at fixed scaling variable x = Q 2 /s. 

Another kinematic domain that is very important at high-energy is given by the (Balitsky- 
Fadin-Kuraev-Lipatov) BFKL limit [3-6], or QCD Regge limit, whereby at fixed Q 2 3> Aq CD , 
s — ► oo. In the BFKL limit, the BFKL evolution in the leading order (LO) governs log(l/a;) 
evolution (at x — > 0) of inclusive processes. Note that the BFKL evolution in the next-to- 
leading order (NLO) [7-10], unlike the LO BFKL [3-5], partly includes GLAPD evolution with 
the running coupling constant of the LO GLAPD, a s (Q 2 ) = An / (3 \og(Q 2 / Aq CD ). 

Therefore, the BFKL and especially the NLO BFKL [7-10] are anticipated to be important 
tools for exploring the high-energy limit of QCD. In particular, this importance arises since 
the highest eigenvalue, u max , of the BFKL equation [3-6, 9, 10] is related to the intercept of 
the Pomeron, which in turn governs the high-energy asymptotics of the total cross-sections: 
a ~ (s/so) aip ~ 1 = (s/so)^ ma ^, where the Regge parameter so defines the approach to the 
asymptotic regime. The BFKL Pomeron intercept in the LO turns out to be rather large: 
oljp — 1 = = 12 log2 (as/n) ~ 0.54 for as = 0.2; hence, it is very important to analyse 
recently calculated NLO corrections [^, §] to the BFKL. 

One of the striking features of the NLO BFKL analysis |J is that the NLO value for 



the intercept of the BFKL Pomeron, improved by the BLM procedure |11|| , has a very weak 
dependence on the gluon virtuality Q 2 : oljp — 1 = io™lo — 0.13 - 0.18 at Q 2 = 1 - 100 
GeV 2 . This agrees with the conventional Regge theory where one expects universal intercept 
of the Pomeron without any Q 2 -dependence. The minor Q 2 -dependence obtained leads to 
approximate conformal invariance. 

There have recently been a number of papers which analyse the NLO BFKL predictions [12- 
16]. Also, a lot of work should be done to clarify the very important issue of the factorisation 
properties of the BFKL regime [17-23]. 

As a phenomenological application of the NLO BFKL improved by the BLM procedure, 
with its effective resummation of the conformal- violating /5o-terms into the running coupling in 



all orders of the perturbation theory, one can consider the gamma-gamma scattering [24, 10 
This process is attractive because it is theoretically more under control than the hadron-hadron 
and lepton-hadron collisions, where nonperturbative hadronic structure functions are involved. 
In addition, for the gamma-gamma scattering the unitarisation (screening) corrections due to 
multiple Pomeron exchange would be less important than in hadron collisions. 

The gamma-gamma cross sections with the BFKL resummation in the LO were considered 
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Figure 1: Virtual gamma-gamma total cross-section by the BFKL Pomeron versus L3 Col- 
laboration data at energies: (a) 183 GeV and (b) 189 GeV of e + e~ collisions. Solid curves 
correspond to NLO BFKL in BLM; dashed: LO BFKL; and dotted: LO Born contribution. 
Two different curves are for two different choices of the Regge scale: Sq = Q 2 /2 and Sq — 2Q 2 . 



in H, [25|, pq| . In the NLO BFKL case one should obtain a formula analogous to LO BFKL 



In Fig. m we present the comparison of BFKL predictions for LO and NLO BFKL improved 
by the BLM procedure with data [^] from L3 at CERN LEP. The different curves reflect the 
uncertainty of the theoretical predictions with the choice of the Regge scale parameter sq, which 
defines the transition to the asymptotic regime. For the present calculations two variants have 
been chosen s = Q 2 /2 and Sq = 2Q 2 , where Q 2 is the virtuality of the photons. One can 
see from Fig. |l] that the agreement of the NLO BFKL improved by the BLM procedure is 
reasonably good at energies of LEP2 ^Js e + e - = 183 - 189 GeV. One can notice also that the 
sensitivity of the NLO BFKL results to Sq is much smaller than in the case of the LO BFKL. 



It was shown in Refs. 128L 29 that the unitarisation corrections in hadron collisions can 



lead to a value of the (bare) Pomeron intercept higher than the effective intercept value. Since 



the hadronic data fit yields about 1.1 for the effective intercept value ||30| , the bare Pomeron 
intercept value should be above it. Therefore, assuming small unitarisation corrections in 
the gamma-gamma scattering at large Q 2 , one can accommodate the NLO BFKL Pomeron 
intercept value 1.13 - 1.18 || in the BLM optimal scale setting, along with larger unitarisation 



corrections in hadronic scattering p9| , p0| , where they can lead to a smaller effective Pomeron 
intercept value of about 1.1 for hadronic collisions. The above intercept value of the NLO 
BFKL Pomeron is in good agreement with the analysis [|31]] of the diffractive dijet production 
at the Tevatron. 

Another possible application of the BFKL approach can be the collision energy depen- 
dence of the inclusive jet production |32|, |22|. Unlike the case with the selection of most 
forward/backward (Mueller-Navelet) jets [19, 33-38], the usual inclusive jets [20-22, 31, 39, 40] 
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can be more reliable for detectors with the limited acceptance in rapidity. 



The advent of the Fermilab Tevatron and the CERN LHC provides a new testing ground 
for the parton model — the kinematic conditions when the energies of the produced hadrons are 
large enough to be described by perturbation theory and, at the same time, are much smaller 
than the total energy of the collision (BFKL semi-hard kinematics). Because the parton model 
was originally invented and subsequently tested for the hard kinematics, the second condition 
makes it plausible that a substantial modification of the parton model will be needed to describe 
this BFKL semi-hard kinematic region. 

The range of applicability of the QCD-improved parton model is a subject of controversy at 
the moment. There are statements (see, e.g., ||l|, [4^1 ) that the fitting capacity of the conven- 
tional QCD-improved parton model is sufficient to accommodate all the data on deep inelastic 
scattering (DIS) parton structure functions available at small- x kinematics. On the other hand, 



the same data from HERA on DIS structure functions can be described by NLO BFKL |IR [16 



In addition, the data from HERA |^3| and Tevatron |44| on most forward/backward jet pro- 
duction may be interpreted as a manifestation of the BFKL Pomeron |33], |I5fl , which is beyond 
the conventional QCD-improved parton model. The situation is further complicated by the ob- 
servation that the range of applicability of the QCD-improved parton model may be different 
for different observables. In particular, the cross-sections of processes with specific kinematics 
exhibit breakdown of the applicability of finite-order perturbative QCD via the development 
of sensitivity to the choice of the normalisation scale. On the other hand, some dedicated 
combinations (ratios) of cross sections may be less sensitive to the inclusion of the higher-order 
corrections. An example is the scaled cross-section ratio [45-48] since, as follows from Ref. j4l| . 
it is relatively insensitive to the inclusion of the NLO correction. 

Under these circumstances, it is crucial to have qualitative predictions from the conventional 
QCD-improved parton model (without BFKL resummation of the energy logarithms) for the 
new kinematic domain. If the predictions would turn out qualitatively incorrect, a generalisa- 
tion, and a substitute in this kinematical domain of the QCD-improved parton model would 
become indispensable. 



Here we discuss such a prediction, made in Ref. pD |. It is a prediction for the ratio of 



inclusive single jet production at a smaller energy y/sjy of the hadron collision to the one at a 
higher energy ^/sp: 

jr>/ \ s N da ,/s D da 

R(x,y,s N ,s D ) = -j—^\ x iyi s N)/ ~faty( x > y > SD >- W 

Here the cross-section is made dimensionless by the rescaling with the corresponding total 
invariant energy of the collision squared sn(sd). The ratio depends on the (pseudo)rapidity 
y = 1/2 ln(k + / k_) , where k± = E ± /c 3 are the light-cone components of the momentum of the 
produced jet, and on the fraction of the energy x = {k + + kJ)/ ' ^fsl, i = N,D deposited in the 
jet produced (sn is used for the definition of x in the numerator, sd in the denominator, so 
that x varies from zero to unity for both energies). Note that this scaling variable coincides in 
the centre-of-mass system with xr = E/E max = 2E / y/s , the radial Feynman variable, or for 
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Figure 2: Inclusive single-jet production: scaled cross-section ratio at the Tevatron. BFKL f32| , 
|22]j and GLAPD |49| predictions versus Tevatron data: (a) CDF data [46], 
data from CDF and D0 Collaborations IBS . 



(b) preliminary 



y — 0, i.e. for (6cms — 7r /2), the scaling variable becomes the transverse Feynman variable: 
x = x± = 2E±/y/s. 



The ratio R (taken at y ~ 0, ie. for jets perpendicular to the collision axes) was used in 



Refs. 45, 46, fi7L [48] as a means to test QCD predictions for scaling violations: 



R(x,y = 0, s N ,s D ) = 




Note that without scaling violations the ratio R is exactly unity. At fixed y and x± the 
dependence of R from sn, sd comes from the presence of the fundamental QCD scale, Aq C d, 
in the running coupling and in the parton distribution functions. 

To be more explicit, we remark that the QCD-improved parton model, based on the fac- 
torisation theorem for hard processes and the GLAPD log Q 2 -e volution, presents the inclusive 
jet scaled cross-section in hadron collsions as 



E 4 ± ^(x ± ,y = 0,s)= f 1 f 1 dx A dx B F A (x A ,Q 2 )F B (x B ,Q 2 )Ei-^5(s + i + u), (4) 

P JXA,min Jx B ,miri 7T CLl 

where s, t and u are the Mandelstam variables for the partonic subprocess, the scale of the hard 
partonic subprocess — t = Q 2 ~ E\ ~ x\s, F A and F B are parton distribution functions with 
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the GLAPD evolution following from perturbative as{Q 2 ) expansion, and the scaled partonic 
subprocess is 



E 1-JI ~ a|(Q 2 )[l + C NLO a s + ...] = a 2 s (x ± s 2 )[l + C NLO a s + ...]. 
at 

Hence, within the QCD-improved parton model, the scaled cross-section ratio for inclusive 
jet production at fixed x and y is the dimensionless function of a.g- The GLAPD scaling 
violation due to the interacting QCD partons appears as the logarithmic^ dependence on the 
total energy of collision through the coupling constant as- We note here that the BFKL leads 
to a power-like scaling violation, the strength of which depends on the Regge scale sq. 



Perturbative QCD calculations of Ref. with hard kinematics (Q ~ s) predict for R 
at y = a steep increase around the value of 1.8 - 1.9 for x growing in the range above 0.1 
(for the case ^/sN/^/s^J = 0.63 TeV/1.8 TeV ). For moderate x, the prediction is in reasonable 
agreement with CDF data |46|, f47|. For x < 0.1, calculations are above the preliminary data of 
CDF |H7|. This was one of the reasons for the conclusion of Ref. [K9T1 that NLO GLAPD 153 



with hard kinematics is insufficient to describe the absolute cross section of jets with transverse 
energy less than 50 GeV within an accuracy of 10%. It was shown in Ref. that resummation 
of the energy logarithms, i.e. BFKL, restores the agreement between theory and experiment. 



In Ref. [5(J we have found the following result: the QCD-improved parton model predicts 



that R is not a monotonic function of its arguments, i.e. the single-jet production cross-section, 
if measured in the natural units of the same cross section taken at another (higher) energy of 
the collision, has extrema. Namely, it has minima ("dips"): there is a value of x for each y with 
the smallest ratio of jets produced. The reason this fact was overlooked is that for y = (the 
only value for which the calculations were reported earlier) the minimum is at a value of x too 
small to be inside the acceptance of the existing detectors (xdi P (y = 0) < 0.01 at the Tevatron). 

Fig. |3]a presents the ratio for energies 0.63 TeV/1.8 TeV at the Fermilab Tevatron, Fig. |3|b 
for energies 6 TeV/14 TeV at the CERN LHC. Each curve on the plots presents the dependence 
of the ratio on x at different values of the rapidity y. Each curve ends at a lowest value of x 
where as(Q 2 ) has the value of about 0.5 (it corresponds to Q = 0.7 GeV, and Q was taken to 
be half of the transverse energy of the jet produced). For lower values of x, perturbative theory 
becomes unreliable because the coupling approaches unity. 

There is another lower bound on the values of x at which our plots make sense, because 
there is a lowest energy for which the jet may be resolved. This energy is accepted now to be 
around 5 GeVQ, which corresponds to x > 0.016 for the Tevatron (Fig. [|a), and to x > 0.0017 
for the LHC (Fig. |b). 

1 Indeed, the small- a; asymptotics of the GLAPD evolution gives a growth of parton structure functions that 
is faster than any power of a logarithm, but slower than any power — so-called double-logarithmic asymptotics 



m h. Buy- 

2 At the HERA lepton-hadron collider jets are resolved from E± = 3 GeV, and at the Tevatron hadron- 



hadron collider the CDF Collaboration tags jets from E± = 8 GeV p4| 



5 




Figure 3: Ratio of scaled cross-sections of inclusive single-jet production: (a) for Fermilab 
energies; (b) for the LHC energies. 



There is an important issue concerning the accuracy of the present LO calculation. The most 
important advantage of the scaled cross-section ratio is that this is the ratio of two perturbative 
series with the same coefficients and with different scales in the running coupling. These scales 
are defined by the two initial collision energies at fixed scaling variable. One can show that the 
theoretical accuracy of the ratio in LO of perturbative QCD is not less than the accuracy of 
NLO calculations for absolute cross-sections. 

The minima in Fig. [| originate from a competition between the running of the parton 
distribution functions and the running of the coupling constant. Namely, the ratio with frosen 
parton distribution functions is decreasing monotonously (this tendency is realised at small x), 
while the one with frosen coupling constant is growing monotonously (which is realised for x 
larger than the position of the minimum). 

We suggest the following potential implications of the minima we have predicted with the 
parton model: (i) If one observes the minima experimentally, one employs the orthodox QCD- 
improved parton model and tries to account for observed positions and depths of the minima 
by taking into account higher-order corrections, in particular, resummation of the energy log- 
arithms, (ii) If one does not observe the minima experimentally, more radical changes are 
motivated, such as an alternative model of the effective constituents inside the hadrons for the 
BFKL semi-hard asymptotics. One example might be the colour dipole model 155 . 



Finally, we comment on the possibility of searching for the minima at the Fermilab Tevatron 
and at the CERN LHC: positions of the minima for the Tevatron energies (Fig. |3]a) seem to 
be reached by both D0 and CDF detectors. The minima of the LHC plot (Fig. |||b) seem to be 
well inside the acceptance of, for example, the FELIX [56], the ALICE ]5J] and the CMS [58 
detectors. 

We take the ratio of 6 TeV/14 TeV for the LHC, because, in addition to 14 TeV pp collisions, 
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lead-lead collisions at the LHC are planned with the collision energy of 6 TeV per nucleon- 
nucleon collision. Since nuclear collisions bring in nuclear effects, which can distort our predicted 
curves, we also considered the ratio 6 TeV/100 TeV (see Ref. 



Further consideration should be given to deciding which pair of energies and value of rapidity 
are most convenient for an experimental search of the minima. Also, more work is needed 
to make quantitative predictions on the locations and the shapes of the dips with the NLO 
corrections taken into account. It is interesting to observe a resemblance of the dips presented 
here with the nonmonotonic behaviour of parton distribution functions in DIS ||59|| . 

It is also worth noting that, in the case of nuclear collisions, the effects of initial nuclear par- 
ton distributions (small- x EMC-effect [pQ| ) and dynamical effects such as quark-gluon plasma, 
jet quenching, etc. pT| , phase structure of the QCD vacuum [B2^ |63| will demand special 
consideration. 

Before reaching conclusions, we would like to note that many of the above ideas can be 
studied also for the case of heavy-quarkonium production, where similar phenomena should be 



present 64 



To sum up, we find a new qualitative prediction of the QCD-improved parton model for 
hadron collisions and suggest its use to test the applicability of the parton model for certain 
regions of high-energy hadron collisions. Study of the scaling violation of the scaled cross 
section ratio can reveal such new dynamic effects as the BFKL asymptotics. 

We thank G. Altarelli, V. P. Andreev, J. Ellis, A. De Roeck, V. S. Fadin, J. H. Field, 
M. Kienzle-Focacci, C.-H. Lin, L. N. Lipatov, J.-W. Qiu, V. A. Schegelsky, A. A. Vorobyov 
and M. Wadhwa for useful discussions, and, also, S. Vascotto for reading of the manuscript. 
VTK thanks the Organizing Committee of the Xth Quantum Field Theory and High Energy 
Physics Workshop and the CERN Theory Division for their warm hospitality. This work was 
supported in part by the Russian Foundation for Basic Research, Grant No. 00-02-17432, 
the NATO Science Programme, Collaborative Linkage Grant No. PST.CLG. 976521, and the 
U.S. Department of Energy, Contract No. DE-FG02-87ER40371, Division of High Energy and 
Nuclear Physics. 



References 

[1] D. Amati, R. Petronzio and G. Veneziano, Nucl. Phys. B140, 54 (1978); ibid. B146, 
29 (1978); S. B. Libby and G. Sterman, Phys. Rev. D18, 3252 (1978); A. H. Mueller, 
Phys. Rev. D18, 3705 (1978); A. V. Efremov and A. V. Radyushkin, JINR-E2-11725, 
11726, 11849, Dubna (1978), Teor. Mat. Fiz. 44, 17, 157, 327 (1980) [J. Theor. Math. 
Phys. 44, 573, 664, 774 (1981)]; R. K. Ellis, H. Georgi, M. Machachek, H. D. Politzer and 
G. G. Ross, Nucl. Phys. B152, 285 (1979); J. C. Collins, D. E. Soper and G. Sterman, 
Nucl. Phys. B261, 104 (1985); G. T. Bodwin, Phys. Rev. D31, 2616 (1985); ibid. (E) D34, 
3932 (1986); for a review, see J. C. Collins, D. E. Soper and G. Sterman, in: Perturbative 
QCD, ed. A. H. Mueller (World Scientific, Singapore, 1989) p. 1. 



7 



V. N. Gribov, L. N. Lipatov, Yad. Fiz. 15, 781, 1218 (1972) [Sov. J. Nucl. Phys. 15, 438, 
675 (1972)]; L. N. Lipatov, Yad. Fiz. 20, 181 (1974) [Sov. J. Nucl. Phys. 20, 94 (1975)]; 
G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977); Yu. L. Dokshitzer, ZhETF 73, 
1216 (1977) [Sov. Phys. JETP 46, 641 (1977). 

E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Phys. Lett. 60B, 50 (1975); ZhETF 71, 
840 (1976) [Sov. Phys. JETP 44, 443 (1976)]; ibid. 72, 377 (1977) [45, 199 (1977)]; L. N. 
Lipatov, Yad. Fiz. 23, 642 (1976) [Sov. J. Nucl. Phys. 23, 338 (1976)]; L. N. Lipatov, 
ZhETF 90, 1536 (1986) [Sov. Phys. JETP 63, 904 (1986)]. 

I. I. Balitsky and L. N. Lipatov, Yad. Fiz. 28, 1597 (1978) [Sov. J. Nucl. Phys. 28, 822 
(1978)]. 

L. N. Lipatov, Phys. Rep. C286, 131 (1997); in: Perturbatwe QCD, ed: A. H. Mueller 
(World Scientific, Singapore, 1989) p. 411. 

L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rep. C100, 1 (1983). 

V. S. Fadin and L. N. Lipatov, Phys. Lett. 429B, 127 (1998); G. Camici and M. Ciafaloni, 
Phys. Lett. 430B, 349 (1998); and references therein. 

A. V. Kotikov and L. N. Lipatov, DESY-00-059 (2000), |hep-ph/ 0004008 . 

S. J. Brodsky, V. S. Fadin, V. T. Kim, L. N. Lipatov and G. B. Pivovarov, 
Pisma ZhETF 70, 161 (1999) [JETP Lett. 70, 155 (1999)], |Eep^ph/990 122S . 

V. T. Kim, L. N. Lipatov and G. B. Pivovarov, IITAP-99-013 and IITAP-99-014, Ames 
(1999), |hep-ph/9911228| and |hep-ph/9911242 . 



S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D28, 228 (1983). 

V. S. Fadin and L. N. Lipatov, Proc. Theory Institute on Deep Inelastic Diffraction, 
ANL, Argonne, September 14 - 16, 1998; C. R. Schmidt, Phys. Rev. D60, 074003 (1999); 
J. R. Forshaw, D. A. Ross and A. Sabio Vera, Phys. Lett. 455B, 273 (1999). 

G. Salam, JHEP 9807, 019 (1998). 

M. Ciafaloni, D. Colferai and G. P. Salam, Phys. Rev. D60, 114036 (1999); M. Ciafaloni 
and D. Colferai, Phys. Lett. 452B, 372 (1999). 

R. S. Thorne, Phys. Rev. D60, 054031 (1999). 

G. Altarelli, R. D. Ball and S. Forte, CERN-TH/2000-010 (2000), |hep-ph/000lT57| , pre- 
sented at 17th Autumn School "QCD: Perturbative or Nonperturbative?" (AUTUMN99), 
Lisbon, Portugal, September 29 - October 4, 1999; CERN-TH/99-317 (1999), |hepl 
ph/9911273 . 



[17] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B366, 135 (1991); J. C. Collins 
and R. K. Ellis, Nucl. Phys. B360, 3 (1991); E. M. Levin, M. G. Ryskin, Yu. M. Shabelsky 
and A. G. Shuvaev, Yad. Fiz. 54 1420 (1991) [Sov. J. Nucl. Phys. 54, 867 (1991)]. 



8 



[18] I. Balitsky, Nucl. Phys. B463, 99 (1996); Phys. Rev. Lett. 81, 2024 (1998); Phys. Rev. 
D60, 014020 (1999). 

[19] A. H. Mueller and H. Navelet, Nucl. Phys. B282, 727 (1987). 

[20] V. T. Kim and G. B. Pivovarov, Phys. Rev. D53, 6 (1996); in: Particle Theory and 
Phenomenology, eds. K. Lassila et al. (World Scientific, Singapore, 1995), p. 366. 

[21] V. T. Kim and G. B. Pivovarov, Phys. Rev. D54, 725 (1996). 

[22] V. T. Kim and G. B. Pivovarov, Proc. of the Vllth Blois Workshop on Elastic and Diffrac- 
tive Scattering: Recent Advances in Hadron Physics (EDS97), Seoul, Korea, June 10 - 14, 
1997, eds: K. Kang, S. K. Kim and C. Lee (World Scientific, Singapore, 1998) p. 236. 

[23] V. T. Kim and G. B. Pivovarov, Phys. Rev. Lett. 79, 809 (1997). 

[24] S. J. Brodsky, V. S. Fadin, V. T. Kim, L. N. Lipatov and G. B. Pivovarov, work in 
preparation. 

[25] S. J. Brodsky, F. Hautmann and D. E. Soper, Phys. Rev. D56, 6957 (1997); 
Phys. Rev. Lett. 78, 803 (1997); ibid. (E) 79, 3544 (1997). 

[26] J. Bartels, A. De Roeck and H. Lotter, Phys. Lett. 389B, 742 (1996); A. Bialas, W. Czyz 
and W. Florkowski, Eur. Phys. J. C2, 683 (1998); M. Boonekamp, A. De Roeck, C. Royon 
and S. Wallon, Nucl. Phys. B555, 540 (1999); J. Kwiecihski and L. Motyka, Acta Phys. 
Pol. B30, 1817 (1999). 

[27] L3 Collaboration, M. Acciari et al, Phys. Lett. 453B, 333 (1999); L3 Collaboration, P. 
Achard, presented at Int. Conf. on the Structure and Interactions of the Photon (Pho- 
ton'99), Freiburg, Germany, May 23 -27, 1999, |hep-ex/99070T6| ; L3 Collaboration, M. 
Kienzle-Focacci, presented at the VHIth Blois Workshop on Elastic and Diffractive Scat- 
tering (EDS99), Protvino, Russia, June 27 - July 2, 1999; M. Wadhwa, presented at the 
19th Int. Conf. on Physics in Collision (PIC'99), Ann Arbor, MI, USA, June 24 - 27, 1999, 
hep-ex/9909001 . 



[28] M. S. Dubovikov, B. Z. Kopeliovich, L. I. Lapidus and K. A. Ter-Martirosian, Nucl. Phys. 
B123, 147 (1977); B. Z. Kopeliovich and L. I. Lapidus, ZhETF 71, 61 (1976) [Sov. Phys. 
JETP 44, 31 (1976)]. 

[29] A. B. Kaidalov, L. A. Ponomarev and K. A. Ter-Martirosyan, Yad. Fiz. 44, 722 (1986) 
[Sov. J. Nucl. Phys. 44, 468 (1986)]. 

[30] J. R. CudeU, K. Kang and S. K. Kim, Phys. Lett. 395B, 311 (1997); 

J. R. CudeU, A. Donnachie and P. V. Landshoff, Phys. Lett. 448B, 281 (1999). 

[31] B. Cox, J. Forshaw and L. Ldnnblad, JHEP 10, 023 (1999). 

[32] V. T. Kim and G. B. Pivovarov, Phys. Rev. D57, 1341 (1998). 



9 



[33 



V. Del Duca and C.R. Schmidt, Phys. Rev. D49, 4510 (1994); 
W. J. Stirling, Nucl. Phys. B423, 56 (1994). 



J. Bartels, V. Del Duca, A. De Roeck, D. Graudenz and M. Wiisthoff, Phys. Lett. 384B, 
300 (1996). 

K. J. Eskola, A. V. Leonidov and P. V. Ruuskanen, Nucl.Phys. B481, 704 (1996). 

C. R. Schmidt, Phys. Rev. Lett. 78, 4531 (1997); 

L. H. Orr and W. J. Stirling, Phys. Rev. D56, 5875 (1997). 

D0 Collaboration, A. Abachi et ai, Phys. Rev. Lett. 77, 595 (1996). 

D0 Collaboration, B. Abbott et ai, the Int. Europhysics Conf. on High Energy Physics, 
Jerusalem, Israel, August 19 - 26, 1997. 

M. Braun and D. Treleani, Eur. Phys. J. C4, 685 (1998). 



R. Peschanski and C. Royon, hcp-ph/0003195 ; J. G. Contreras, R. Peschanski and C 



Royon, DAPNIA-SPP-00-04, Saclay (2000), [hep-ph/ 0002057 



R. D. Ball and S. Forte, Phys. Lett. 335B, 77 (1994); ibid. 405B, 317 (1997). 

F. Barreiro, C. Lopez and F. J. Yndurain, Z. Phys. C72, 561 (1996); F. J. Yndurain, 
FTUAM 96-19, Madrid (1996), |hep-ph/9605265| . 

ZEUS Collaboration, J. Breitweg et al, Eur. Phys. J. C6, 239 (1999); 
HI Collaboration, C. Adloff et al, Nucl. Phys. B538, 3 (1999). 



D0 Collaboration, B. Abbott et al, FERMILAB-PUB-99-363-E, Batavia (1999), [Eepj 
ex/9912032. 



UA2 Collaboration, J.A. Appel et al, Phys. Lett. 160B, 349 (1985); 
UA1 Collaboration, G. Arnison et al, Phys. Lett. 172B, 461 (1986); 
CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 62, 613 (1989). 

CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 70, 1376 (1993). 

CDF Collaboration, T. Devlin, Proc. of the XXVIII ICHEP'96, Warsaw, Poland, July 25 
- 31, 1996, eds. Z. Ajduk and A. K. Wroblewski (World Scientific, Singapore, 1997); A. A. 
Bhatti, Proc. of the DPF'96 Meeting, Minneapolis, USA, August 10 - 15, 1996 (World 
Scientific, Singapore, 1997). 

[48] J. Huston, Proc of the XXIX ICHEP'98, Vancouver, Canada, July 23 - 29, 1998, eds. A. 
Astbury, D. Axen and J. Robinson (World Scientific, Singapore, 1999) Vol. 1, pp. 283 - 
302, |hep-ph/990135l ; 

CDF Collaboration (preliminary), see: A. Akopian, PhD thesis, Rockefeller Univ. (1999), 
fittp : / /www — cdf . f nal.gov / physics / new / qcd / QC D .html\ ; 



D0 Collaboration (preliminary), see: J. Krane, PhD thesis, Nebraska Univ. (1999), 
fittp : I /www — dO.fnal.gov/results/publicationsJalks/thesis/thesis.htrnl] 



10 



[49] S. D. Ellis, Proc. of the XXVIIIth Rencontres de Moriond, Les Arcs, France, March 20 - 
27, 1993, ed. J. Tran Thanh Van (Editions Frontieres, 1993) Vol. 2, p. 235. 

[50] V. T. Kim, G. B. Pivovarov and J. P. Vary, Foundations of Physics 30, 519 (2000), [hep] 
| ph/9709301 . 

[51] A. De Riijula, S. L. Glashow, H. D. Politzer, S. B. Treiman, F. Wilczek and A. Zee, Phys. 
Rev. D10, 1649 (1974). 

[52] F. Olness and W.-K. Tung, Int. J. Mod. Phys. A2, 1413 (1987); W.-K. Tung, Nucl. Phys. 
B315, 378 (1989). 

[53] S. D. Ellis, Z. Kunszt and D. E. Soper, Phys. Rev. Lett. 64, 2121 (1990); 

F. Aversa, M. Greco, P. Chiappetta and J.Ph. Guillet, Phys. Rev. Lett. 65, 401 (1990); 
W. Giele, E. W. N. Glover and D. A. Kosower, Phys. Rev. Lett. 73, 2019 (1994). 

[54] CDF Collaboration, F. Abe et ai, Phys. Rev. D57, 67 (1998). 

[55] A. H. Mueller, Nucl. Phys. B415, 373 (1994); 

N. N. Nikolaev and B. G. Zakharov, Phys. Lett. 327B, 149 (1994). 

[56] K. Eggert and C. Taylor, FELIX: a Full Acceptance Detector for the CERN LHC, 
CERN-PPE-96-136 (1996); FELIX Collaboration, E. Lippmaa et al, Letter of Intent, 
CERN/LHCC 97-45 (1997). 

[57] ALICE Collaboration, Technical Proposal, J. Schukraft et ai, CERN/LHCC 95-71, Geneva 
(1995). 

[58] CMS Collaboration, Technical Proposal, M. Delia Negra et al, CERN/LHCC 94-38, 
Geneva (1994). 

[59] P. Desgrolard, L. Jenkovszky and F. Paccanoni, Eur. Phys. J. C7, 263 (1999); P. Desgro- 
lard, L. Jenkovszky, A. Lengyel and F. Paccanoni, Phys. Lett. 459B, 265 (1999). 

[60] EMC Collaboration, J. J. Aubert et ai, Nucl. Phys. B293, 740 (1987); 
M. Arneodo, Phys. Rept. C240, 301 (1994). 

[61] X. N. Wang, Phys. Rep. C280, 287 (1997). 

[62] G. 't Hooft, Nucl. Phys. B153, 141 (1979); M. Luscher, Nucl. Phys. B219, 233 (1983); P. 
van Baal, Nucl. Phys. B369, 259 (1992). 

[63] V. T. Kim, V. A. Matveev, G. B. Pivovarov and J. P. Vary, IITAP-99-012, Ames (1999), 
|hep-ph/9910T57 . 

[64] V. T. Kim, G. B. Pivovarov, H. S. Song and J. P. Vary, work in progress. 



11 



